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ABSTRACT 



We present the results of deep I-band imaging of two BL Lacerate objects, RGB 0136+391 and PKS 0735+178, during an epoch 
when the optical nucleus was in a faint state in both targets. In PKS 0735+178 we find a significant excess over a point source, which, 
if fitted by the de Vaucouleurs model, corresponds to a galaxy with 1= 18.64 ±0.11 and r e jf = 1.8 ±0.4 arcsec. Interpreting this galaxy 
as the host galaxy of PKS 0735+178 we derive z = 0.45 ± 0.06 using the host galaxy as a "standard candle". We also discuss the 
immediate optical environment of PKS 0735+178 and the identity of the Mgll absorber at z = 0.424. Despite of the optimally chosen 
epoch and deep imaging we find the surface brightness profile of RGB 0136+391 to be consistent with a point source. By determining 
a lower limit for the host galaxy brightness by simulations, we derive z > 0.40 for this target. 

Key words. Galaxies:active - BL Lacertae objects: individual: RGB 0136+391, PKS 0735+178 - Galaxies: nuclei 



1. Introduction 

The optical spectrum of BL Lacertae objects (BL Lacs) is domi- 
nated by a featureless synchrotron continuum, which makes the 
determination of their redshift by spectroscopy very challeng- 
ing. Redshift determination of quasars and BL Lacs is usually 
based on emission lines originating in the narrow and broad line 
regions and/or on the absorption features of the host galaxy. 
However, in BL Lacs the emission lines are intrinsically weak 
and like the host galaxy absorption lines they often become so 
weak relative to the continuum that they fall below the noise 
level and become unrecognizable. 

It is therefore not surprising that many BL Lac samples still 
have a high fraction of targets without spectroscopically deter- 
mined redshifts. For instance, in one of the most studied sam- 
ples, the 1 Jy sample (IStickel et al.ll 199 II) . 23 of 37 targets have a 
spectroscopically secured redshift with the remaining being un- 
known, lower limits from an absor ption line system or un certain 
due to a detection of a single line (Rector & Stocke 2001). Large 
(8 m and up) telescopes ar e often needed to secure the redshift 
(e.g. ISbarufatti et al.ll2009l) . but even then there often remains a 
fraction of targets with unknown redshifts. 

ISbarufatti et al.l d2005l) discussed in detail the method of es- 
timating the redshift of BL Lacs through detection of their host 
galaxies. They showed that the absolute magnitude distribution 
of BL Lac host galaxies is roughly Gaussian with an average 
of Afjj° st = -22.8 and <x = 0.5 mag. This opens the possibil- 
ity of estimating the redshift using the host galaxy as a "stan- 



* Based on observations made with the Nordic Optical Telescope, op- 
erated on the island of La Palma jointly by Denmark, Finland, Iceland, 
Norway, and Sweden, in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica de Canarias. 



dard candle" if sufficiently deep optical or NIR images can be 
obtained. This method has been applied to several BL Lacs to 
estim ate their redshifts (Nils son et al.ll2008l : iMeisner & Ro mani 
l2010h . Given the width of the host galaxy luminosity distribu- 
tion, the imaging redshifts are not as accurate as spectroscopic 
redshifts. The accuracy of the method is +0.05 in redshift, if high 
signal-to-noise data are available, but in practice the redshift er- 
ror is usually larger du e to noise in derived host galaxy properties 
(ISbarufatti et al.ll2005l) . 

In this paper we report the results of deep optical imaging of 
two y-ray bright BL Lacs, RGB 0136+391 and PKS 0735+178, 
which were detected b y the Fermi-LAT already during first three 
months of operations dAbdo et al]|2009l) . For y-ray emitting BL 
Lac objects the redshift estimation is particularly important and 
there are large dedicated spectroscopy programs for determin- 
ing redshift for Fermi-LAT detected AGN. Still 56% of the 
LAT detected BL Lac objects have no redshift measurement 
(lAckermann et al.ll201 ll) . which makes population studies diffi- 
cult . RGB 0136+391 is particularly interesting, because of its 
very hard y-ray spectrum (spectral index 1.69) in the Fermi- 
LAT range and the very recent discovery of VHE y-rays from 
the source (M azinll2012b . VHE y-rays do not travel unaffected 
in the intergalactic space, but are absorbed by the extragalac- 
tic background light (EBL). The EBL leaves an imprint on the 
observed energy spectrum and if the redshift of the emitting ob- 
ject is known, the observed energy spectrum can be used to set 
constrains to density and evolution of EBL with redshift (e.g . 
Stecker et al lll996t [Aharonian et al.ll2006HMazin & Raue 2007; 



Mever et al 1l2012l) . Alternatively, if the redshift is not known the 
observed VHE y-ray spectrum combined with the GeV range 
spectra can also be us ed to derive independent es timation for the 
redshift of the object dPrandini et al .11201 Oi 1201 lh . 
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RGB 0136+391 was f i rst ide ntified as a BL Lac object by 
lLaurent-M uehleise n et al. I (Il998h who saw a fe atureless opti- 
cal spectrum. Following attempts (|Wei et al.lfl99 9; Bauer etaU 
2000; ICaccianiga etai1l2002t iPiranomonte et al.ll2007l) have all 



failed t o detect a ny spectral featu res and the redshift remains un- 
known. Piranomonte et al. ( 2007 ) derived a lower limit z > 0. 1 93 
based on the absence of host ga laxy features in the s pectrum. In 
the only imaging attempt so far. lNilsson et al.l(l2003l) did not de- 
tect the host galaxy in a 1200s R-band exposure. 

PKS 0735+178 w as classified and a BL Lac object by 
ICarswell et al.J (1 19741) . who detected no emission lines, but a 
strong /12798A Mgll absorption system at z=0.42 4 providing a 
lower limit for the redshift. Later spectroscopy (e . g . IStickel et al.l 
11993b iFalomo et al.ll 1994 iRector & Sto cke 200 1 ) has confirmed 
the absorption system, but no emission features or host galaxy 
absorption lines have bee n detected and the redshift remains un- 
known. Optical imaging (Hutchings et al. 1988; Abraham et al 
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Stickeletal. 1993; Falomo & Ulrichl 120001: IScarpa et al . 
Pursimo et al 2002) has also failed to detect the host 



host galaxy based on a 440 s HST image through the F702W 
filter. 

Despite of the str ong (equivalent width W b s ~ 3. 6 A, 
IRector & Stockd l200ll) Mgll absorption line , the absorbing 
galaxy has remained hidden. [Hutchings et al. (1988) reported 
a bridge of R-band emission from PKS 0735+178 towards the 
galaxy 7 arc sec NW with a SB o f 26.0-25.2 mag/sq. arcsec. Later 
imaging bv IStickel etaII(fT99l did not see this bridge and the 
galaxy 7 arcsec NW was measured to have z = 0. 645. There are 
two p ossible detections of the absorbing galaxy: iPursimo et al.l 
( 1999) reported a faint (R = 21 .9) elongated s tructur e 3 arcsec 
NE of PKS 0735+178 and IFalomo & Ulrichl (l2000h saw faint 
emission 3.5 arcsec E of 0735+178 suggesting it could be the 
intervening z=0.424 system, although they couldn't rule out the 
possibility that it is just a background target. 

The brightness of the BL Lac nucleus relative to the host 
galaxy often makes it very difficult to determine the host galaxy 
parameters (magnitude and effective radius) accurately or to de- 
tect the host galaxy in the first place. Given that the nucleus 
is often highly variable, it is advantageous to wait for the mo- 
ment when it is in a low state. We are currently performing op- 
tical monitoring of -60 blaza rs, including R GB 0136+391 and 
PKS 0735+178 (seeQ and e.g. lAleksic etaDl201 ll for results of 
this monitoring). We are running a parallel program at the Nordic 
Optical Telescope (NOT) to perform deep I-band imaging of a 
sample of gamma-ray emitting blazars or potential gamma-ray 
emitters whenever our monitoring shows that the optical bright- 
ness has fallen below a pre-defined threshold. Instead of observ- 
ing the whole sample at once we wait for the optimal moment. 
To improve the chances of detecting the host galaxy the thresh- 
old level is chosen to represent a very low optical level of each 
target considering the historical light curve. 

In December 2008 the R-band brightness of RGB 0136+391 
was observed to be close to the lowest level since the beginning 
of our monitoring (fall 2002) and in January 201 1 the brightness 
of PKS 0735+178 reached R = 16.7, a ve ry low level consid - 
ering the long-term historical light curve dCiprini et al.ll2007l) . 
Observations were promptly triggered and performed at the NOT 
and the results are reported in this paper. 

Throughout this paper we use the cosmology Hq = 70 km 
s _1 Mpc~', Q.m — 0.3 and £2a = 0.7. 



*6 



/ 



/■ 



■ • 

1* . 



5. 
♦ 



Fig. 1. Summed image of RGB 0136+391 with the stars dis- 
cussed in the text labeled. Field size is 5.7 x 5.8 arcmin. North 
us up and east is to the left. 

2. Observations and data reduction 

RGB 0136+391 was observed at the Nordic Optical Telescope 
(NOT) through an I-band filter with almost uniform transmission 
between 725 and 825 nm on December 7th, 2008. We used the 
ALFOSC camera with a gain of 0.726 e /ADV and readout noise 
of 4.2 e~. The field of view of the 2048 2 E2V chip was ~ 6.5x6.5 
arcmin (0. 1 89 arcsec/pixel). Altogether 131 images with individ- 
ual exposure times between 25 - 100s were obtained resulting in 
a total exposure time of 3830s. In order to make a fringe cor- 
rection image the telescope was moved between each exposure. 
The images were bias-subtracted and flat-fielded with twilight 
flats and a fringe correction image was constructed from the ex- 
posures and subtracted. The images were then registered using 
15 unsaturated stars in the field and coadded. The summed image 
shown in Fig.[T]has a FWHM of 0.61 arcsec. Due to the presence 
of a quarter Moon in the sky the background brightness was 19.2 
mag/sq. arcsec, ~ 0.8 mag brighter than average dark sky on La 
Palmdi The field was calibrated using I-band observations made 
with the 35cm KVA (Kungliga Vetenskapsakademien) telescope 
on La Palma in photometric cond itions on Jan 1 1th, 2012. The 
standard star field PG 0231+055 dLandoltll 19921) was observed 
at the same airmass as the RGB 0136+391 field and the I-band 
magnitudes of stars 1-6 were derived from these observations. 
The same stars were then used to determine the zero point of our 
summed image. The error of the derived zero point is 0.02 mag, 
excluding color effects (see below). 

PKS 0735+178 was observed with the same instrumental 
setup at the NOT on January 18th, 2011. Altogether 55 images 
with exposure times of 150 - 300 s were obtained. The images 
were processed in a way identical to RGB 0136+391. The fi- 
nal coadded image shown in Fig. [2] has a total exposure time 



http://users.utu. fi/~kani/ 1 m/index. html 



2 La Palma technical note 115, 

http://www.ing.iac.es/Astronomy/observing/conditions/skybr/skybr.html 
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Fig. 2. Summed image of PKS 0735+178 with the stars dis- 
cussed in the text labeled. Field size is 5.2 x 5.4 arcmin. North 
us up and east is to the left. 

of 10250s and a FWHM of 0.60 arcsec. The observations were 
influenced by almost full (~ 90%) Moon 21 degrees away from 
PKS 0735 + 178. As a result the sky background is relatively high 
(17.7 mag/sq. arcsec) and there is scattered light present in the 
raw images resulting in an uneven background. The fringe re- 
moval procedure mostly eliminates the scattered light but leaves 
a small (0. 1 % peak-to-peak amplitude) undulation of the back- 
ground at large spatial scales, barely visible in Fig. [2] Therefore, 
the background of PKS 0735+178 and its PSF stars was sub- 
tracted by fitting a two-dimensional linear slope to the sky pix- 
els surrounding the targets, effectively removing any remaining 
tilt in the backgr ound. Calibration of the field was done through 
star 5, for which lCiprini et al.l J2007I) give I = 15.12 + 0.06. This 
star was unsaturated in several exposures, which were used to 
derive I-band magnitudes for stars 1-3, which were then used to 
determine the zero point of the coadded image. 

Since the transmission curve of the filter use d here closely 
matches the I-band transmission (Bessell 1990), the color ef- 
fects are expected to be small. According to the NOT zero 
point monitoring pagefl the (V-I) color term is -0.03. Two of 
the stars used for calibration, star 4 in RGB 0136+391 and 
star 5 in PKS 073 5+178, have V-I * 0.8 dNilsson "et alJl2007t 
ICiprini eF al. 2007, present work). Assuming the rest of the stars 
used for calibration have similar colors an d that the host galaxy 
V-I = 2.6 (elliptical galaxy color at z = 0.5. lFukugita et al.lll995l) 
we derive the error of host galaxy magnitudes due to filter mis- 
match to be as 0.05 mag. 

3. Analysis 

We searched for evidence of the host galaxy by fitting two- 
dimensional models to the observed image. The detections are 

3 http://www.not.iac.es/instruments/alfosc/zpmon/ 



further tested by simulations, which are also used to derive error 
bars for the derived parameters. 

The fitting pro cedure is the same as em ployed in our pre- 
vious studies (e.g. lNilsson etal]ll999Ll2008h . In short, the two- 
dimensional model consists of two components, the unresolved 
core and a host galaxy whose surface brightness profile follows 
the de Vau couleurs law (J3 = l/n = 0.25, where n is the profile 
index from lSersiq [l968). The model is described by five param- 
eters: location x-y, assumed to be the same for the core and the 
host galaxy, core magnitude 7 core , host galaxy magnitude /host 
and host galaxy effective (half light) radius r e g. Since the host 
galaxies were known from previous imaging to be very faint, the 
ellipticity and position angle were kept at constant value = 0. 
We used a Levenberg-Marquardt loop to find the set of param- 
eters which minimized the x 1 between the data and the model. 
The models were convolved with an empirical point spread func- 
tion (PSF) derived from stars in the vicinity of the target. In the 
case of RGB 0136+391 we used stars 1-4 and in the case of 
PKS 0735+178 stars 1-3. The fitting proceeded in two phases. 
First we fitted a model consisting of an unresolved component 
only, i.e. the model had three free parameters (x, y and / core ). 
After this the position was fixed and the full core + host galaxy 
model was fit with 7 core , 7h os t and r e g as free parameters. 

If the host galaxy is below detection limit these fits usually 
converge towards r e ff — > or r e ff — » oo. In case we found r e ff > 
we tested the significance of the detection by performing fits to 
simulated data. We created a set of ~100 simulated images of 
the target using the fitted core and host galaxy values. The f3 
value for the surface brightness profile of the simulated galaxy 
was drawn from a Gaussian distribution with average equal to 
0.25 and <x = 0.025. The fitting, however, was made with fixed 
/3 = 0.25. Each model was convolved with a PSF consisting of 
a Moffat /3 = 2.5 profile with a FWHM equal to the observed 
one, an ellipticity randomly drawn from an uniform distribution 
between 0.0 and 0.06 and a position angle drawn from an uni- 
form distribution between 0.0 and 180 degrees. To these images 
we added Gaussian noise representing photon and readout noise, 
scaled to the same level as in the observed data and a random off- 
set drawn from a Gaussian distribution with <x = cr sky /30, where 
cr^y is the sky rms scatter, representing the error made in sky 
level estimation. 

For each simulation we also constructed a second PSF, which 
slightly deviated from the PSF used in convolving the model. 
The ellipticity of the second PSF was drawn from an uniform 
distribution with mean equal to the ellipticity of the first PSF 
and maximum deviation of 0.02. Similarly, the position angle 
deviated maximum 10 degrees from the position angle of the 
first PSF. To this second PSF we also added random photon and 
readout noise plus an offset the simulate the sky determination 
error. The second PSF was then used in the simulated fits as the 
assumed PSF. Using a different PSF for the model creation and 
for the fitting enabled us to roughly simulate the PSF variation 
seen in the images and to produce residuals which qualitatively 
resemble the residuals seen in actual images, although the latter 
tend to be a bit more complex than what these simple simula- 
tions can produce. The limits to the distributions above were se- 
lected to match the residuals quantitatively, i.e. the peak-to-peak 
residuals in the simulations were on average equal to the real 
residuals. Similarly, the average^ 2 in the simulations was close 
to the x 1 of the fit to the real data. 

After the simulations we computed the standard deviation of 
the host galaxy magnitudes cr host and the range of effective radii 
which includes 67% of the simulated values. We didn't use stan- 
dard deviation for the effective radius because the distribution 
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Fig. 3. I-band surface brightness profile of RGB 0136+391 
(thick solid line) and four stars in the field. The stars are la- 
beled in Fig.Q] Insert: Residuals after subtracting the PSF from 
RGB 0136+391. Field size is 14x14 arcsec. The error bars in- 
clude contribution from the empirical scatter in each bin and the 
error in sky background determination. 

was asymmetric with a longer tail towards high values. To con- 
sider a host galaxy detected we required <x host < 0.3 mag and 
that the median effective radius of the simulations is > 2cr reff , 
where cr K s is the lower 67% error bar. The more relaxed crite- 
rion for r e ff was used because r e ff is generally more difficult to 
constrain and requiring r e $ > 3cr reff would exclude host galaxies 
whose surface brightness profile clearly exceeds any reasonable 
PSF variability. 

4. Results 

4.1. RGB 0136+391 

Figure [3] shows the surface brightness (SB) profile of 
RGB 0136+391 together with the profiles of four stars in the 
FOV. The profiles of the stars have been scaled to the same to- 
tal flux as RGB 0136+391. The SB profile of RGB 0136+391 
is totally consistent with a point source, with no excess in 
the outer part of the profile to indicate the presence of a host 
galaxy. Examination of the two-dimensional residuals (Fig. [3] 
insert) yields the same conclusion: no obvious signs of the host 
galaxy above PSF residuals can be seen. The model fits con- 
firm this conclusion: the final values returned by the fits are 
= 15.20+0.01, W = 19.46 + 0.32 and r eB = 1 . 5^ j arcsec, 
where the error bars come from the simulations described above. 
This target thus remains unresolved. 

We then derived a lower limit for the redshift of 
RGB 0136+391 using the following procedure: if the distri- 
bution of BL Lac ho st galaxy magnitudes is as shown by 
ISbarufatti et al.1 d2005l) . then 95% of BL Lac host galaxies are 
brighter than M R = -22.0. We thus created 100 simulated BL 
Lac objects at z = 0.25, 0.30, 0.45 with core magnitude 
/core = 15.20 and /host an d r eS corresponding to M^ ost = -22.0 
and Z? e ff = 10 kpc. We first computed the absolute I-band mag- 



Fig. 4. Surface brightness decomposition of PKS 0735+178. In 
addition to PKS 0735+178 and model surface brightnesses, the 
profiles of stars 1-4 are shown. Insert: the residuals after sub- 
tracting the core + host galaxy model. Field size is 14.2 x 14.2 
arcsec. 

nitude Mj using R - I = 0.70 (iFukugita et al.ll 19951) and the ap- 
parent I-band magnitude nti from 

mi - Mi + DM + K!+Aj- E(z) (1) 

where DM is the distance modulu s, Ki is the I-band K- 
correction from Fukugit a et al.1 (U9 95). Aj — 0.15 is the galac- 
tic extinction dSchlegel et alj 1 19981) and E(z) = 0.84 * z is 
the evolution correction, comput ed using the PEGASE code 
dFioc & Rocca-Volmerangelll997l) by assuming galaxy forma- 
tion 1 1 Gyr ago and passive evolution thereafter. 

At each redshift we then fitted the 100 model images using 
the same procedure as in the error simulations and computed the 
standard deviation of /h os t and r eB . The result of these simula- 
tions was that at z = 0.45 (;«/ = 19.52) the host galaxy became 
unresolved using our criteria. Thus at z < 0.40 we would have 
detected the host galaxy with >95% probability and thus assign 
a conservative lower limit z > 0.40 for RGB 0136+391. 

4.2. PKS 0735+178 

Figure g] shows the SB profile of PKS 0735+178 together with 
stars 1-4. The stellar profiles have been again scaled to the 
same total flux as PKS 0735+178. Contrary to RGB 0136+391, 
there is a clear excess over the stellar profiles. Model fits yield 
/ C o re = 16.19 ± 0.01, /host = 18.64 + 0.07 and r eB = 1.8 + 0.4 
arcsec meeting our criteria for successful detection. Adding the 
error from zero point determination and filter mismatch the er- 
ror of /host becomes 0. 1 1 mag. Given the clear excess and the fit 
results we conclude that we have detected the host galaxy or a 
foreground galaxy (see Sect. [5]). 

The error bars from the simulations do not take into account 
possible systematic errors, which are difficult to predict and/or 
quantify. We have tested the sensitivity of our result with respect 
to our particular choice of software by performing the fits with 
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Fig. 5. 47.4 x 47.4 arcsec field around PKS 0735+178 after sub- 
tracting the scaled PSF from PKS 0735+178 and star 4. The im- 
age has been smoothed with a Gaussian kernel with cr = 1.0 
pix. Galaxies Gl and G2 and the rectangular test feature R are 
discussed in the text. H is the host galaxy of PKS 0735+178. 
The dotted ellipse and circle mar k th e location of the feature s 
discussed bv lPursimo eTaTl (1 19991) and lFalomo & Ulrichl (§000), 
respectively. 



the GALFIT software dPeng et al.ll2002l) . We used the same im- 
age, PSF and sky background as before with the only difference 
being the software. The result from this fit was /host = 18.57 and 
r e ff = 2.0 arcsec, i.e. the host galaxy magnitude differs by 7% 
and the effective radius by 1 1 % from our result. The two pro- 
grams thus provide very consistent results and we do not expect 
the systematic errors due to software choice to dominate over 
statistical noise. 

Assuming that the detected galaxy is the host galaxy we 
est imated the redshift o f PKS 0735+178 using the results 
of ISbarufatti et all d2005l) . Using Eq. [TJ with Aj = 0.068 
dSchlegel et al.ll 19981) and M/ = -22.8 we iteratively determined 
the the redshift consistent with the observed I-band magnitude 
to be z = 0.45 + 0.06, where the error bar includes the contribu- 
tion from host magnitude error and the uncertainty of the method 
itself. 

5. Discussion of PKS 0735+178 



Scarp a et"a"D (|2000) derived an upper limit niR > 20.44 for the 
host galaxy, which together with our result implies R-I > 1.7. 
This is redder than the ex pected R-I = 1.2 for ellipticals at 
z=0.5 (Fukug ita et alJ[l995h . possibly an indication of system- 
atic differences between the two studies, or of dust reddening 
in the host galaxy. If significant dust absorption is present in 
PKS 0735+178, our redshift estimate is too high. 

Using the I-band data obtained here and previous R-band 
data we now discuss the identity of the Mgll absorber. Galaxy 
Gl in Fig. [5] can b e directly ruled out d ue to the redsh i ft z = 
0.64 5 measured by IS ticket et al.1 d!993l) . iPursimo et ail d 19991) 
and iFalomo & Ulrichl d2000) detected weak R-band emission 



3'.'0 NE and 3'.'5 E of PKS 0735+178, respectively, and discussed 
the possibility of these features being the Mgll absorber. Neither 
of these emission features are visible in our PSF-subtracted I- 
band image (Fig. |5). Since our images were obtained in bright 
moonlight and thus not optimal for detecting faint diffuse emis- 
sion, we tested our sensitivity with the following procedure. 

Firstly, since Falomo & Ulrichl (2000) do not give a mag- 
nitude for the feature they detected, we cannot estimate if 
it was detectable in our image. We therefore concentrate on 
the elongated featur e reported 370 NE of PKS 0735+178 by 
IPursimo et al.l d 1999b . To roughly simulate this feature we cre- 
ated a rectangular test object marked with R in Fig. [5] Note that 
we did not create this feature at the correct location to better 
compare with the actual image. The feature has a uniform I- 
band surface brightness of 24.3 mag/sq. arcsec and its dimen- 
sions and integrated magnitude 21.9 correspond to th ose of the 
elongated featured reported by Pursi mo et al. (1999). The test 
feature is just visible o ver the background noise and thus we 
should have detected the lPursimo et al.l (1 19991) feature in our im- 
age if its R-I color > 0. Since the range of R-I colors of E to Im 
galaxi es at z = 0.424 is exp ected to be form 1.1 to 0.4, respec- 
tively dFukugita et al.lll995l). we conclude that we should have 
detected the lPursimo et alJd 19991) feature, if its colors are within 
the range of "no rmal" galaxies . Given the non-detection in our 
I-band image the lPursimo et al.l (1 19991) feature is unlike to be the 
absorber. It could be a line emitting cloud whose emission wave- 
length falls into the R-band but outside the I-band. Of the lines 
bluewards of the R-band, A50Q1 [OIII] is the most likely candi- 
date. This line would place the cloud at z = 0. 13-0.44 when 50% 
transmission limits of the R-band filter are considered. 



Since we do no detect the feature discussed by Pursimo et al 



1999) and the detection probability of the IFalomo & Ulrich 



(2000) feature cannot be estimated we turn our attention to 



galaxy G2. This galaxy lies at projected distance of 872 (46 kpc) 
from PKS 0735+178. Given the rest frame equivalent width of 
the Mgll fine in PKS 0735+178, W = 3.6A/(1 + z ab s) = 2.5 k, 
we ask if this ga l axy co uld be the absorber. For comparison we 
use lZibetti etalj d2007l) . who studied 2844 Mgll absorbers with 
W > 0.8A and 0.37 < z < 1 using the SDSS DR4. We first 
note that Wo = 2.5A implies very strong absorption: 93% of 
the systems studied by Zibetti et alj d2007l) have W < 2. 5 A. If 
we concentrate on the relevant section of the parameter space 
in lZibettietail d2007l) . 0.37 < z < 0.55 and W > 1.58A, we 
note that Mgll absorbers can be found up to projected distances 
of ~ 100 kpc from the background source. In this section of the 
parameter space the probability of finding an Mgl l absorber at 
proje cted distance of 46 kpc is ~ 25% (see Fig. 9 in Zibett i et al.l 
2007). Given that the probability of finding an Mg II absorber at 
a given distance decreases with Wo and the Wo of the Mgll line 
in PKS 0735+178 is at the high end of the distribution, the prob- 
ability of G2 being the absorber is probably considerably lower. 
The I-band magnitude of G2 is 20.4 + 0.1 whic h gives R-I = 
0.5 wh en combined with the R-band magnitude in Pursim o et al.l 
dl999l) . This color is consistent with a late type (Scd) galaxy at 
z = 0.424. The absolute magnitude of G2 at z = 0.424 would be 
M R = -20.9 ~ M* + 0.3. In the light of these results G2 can- 
not be completely ruled out as the absorbing galaxy, although 
finding such a system is relatively unlikely. 

As the next option we consider the possibility that the de- 
tected host galaxy is actually a foreground object and not the real 
host galaxy. If this is the case, then the impact parameter must 
be very close to zero since the detected galaxy appears to be well 
centered on the BL Lac nucleus. To check for possible offsets be- 
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tween the core and the host galaxy we performed the host galaxy 
fit with 7 free parameters, i.e. in addition to 7 core , /host and r e ff also 
the x-y positions of the core and the host galaxy were allowed to 
change freely. The core - host galaxy offset was 0.67 pixels (0.13 
arcsec) in this case, indicating that the core is well centered on 
the host galaxy. Such a system is capable of producin g gravita- 
tional lensing phenomena, but R ector & Stockd (120031) did not 
detect any signs of gravitational macrolensing (Einstein rings or 
multiple images) in their 8.46 and 14.94 G Hz maps with ~ 0.2 
arcsec resolution. Furthermore, the result o f IZibetti et al.l d2007l) 
indicate that strong absorbers are associated with galaxies with 
intense star forming, which should be visible in the optical spec- 
trum as narrow emission lines. These findings do not disprove 
the hypothesis that the detected galaxy is an intervening system 
with a small impact parameter, but they make it less attractive 
with respect to the alternative, that we have detected the actual 
host galaxy. 

Lastly we note that the derived redshift 0.45 + 0.06 is 
within errors equal to the absorption redshift 0.424. It could 
thus be possible that the absorption occurs in the host galaxy 
of PKS 0735+178 itself. This would imply the PKS 0735+178 
host galaxy to be a gas-rich system since strong Mg II absorption 
is often associated with dam ped Lyman-alpha systems (DLAs) 
with high HI gas content dPeroux et al.l I2004T) . Further obser- 
vations are clearly needed in order to determine which, if any, 
of the above scenarios is the correct one. A spectrum of G2 is 
needed in order to see if its redshift is compatible with being the 
absorbing galaxy. An ultraviolet spectrum of PKS 0735+178 de- 
tecting the Lyman-alpha line would also shed new light on the 
nature of the absorber by measuring of the HI column density. 

6. Conclusions 

We have presented deep I-band imaging of two BL Lacertae ob- 
jects, RGB 0136+391 and PKS 0735+178, obtained during a 
deep minimum in the optical light curve. Despite of a deep ex- 
posure obtained under excellent seeing conditions (0V6) the host 
galaxy of RGB 0136+391 remains unresolved. We derive z > 
0.40 for RGB 0136+391 using simulated host galaxies at differ- 
ent redshifts with Mr = -22.8 and = 10 kpc. This makes 
RGB 0136+391 one of the most distant VHE y-ray sources de- 
tected up to date. 

The surface brightness profile of PKS 0735+178 reveals an 
excess over a point source, which, if fitted by a de Vaucouleurs 
profile, corresponds to a galaxy with I = 18.64 + 0.1 1 and r e ff = 
1.8 + 0.4 arcsec. Under the assump tion that BL Lac ho st galaxies 
can be used as standard candles (ISbarufatti et al.l2005l) we derive 
z = 0.45 + 0.06 for PKS 0735+178. 

We also discuss the identity of the galaxy responsible for the 
strong Mgll absorption line in PKS 0735+178 at z abs = 0.424 
using previous R-band imaging and the I-band data in this paper. 
We do not detect the two faint objects suggested earlier to be 
responsible for the absorption or any other obvious candidates 
in the vicinity of PKS 0735+178. We identify a galaxy 8'.' 2 of 
PKS 0735+178 (46 kpc projected distance) as a candidate for 
the absorber, although the probability of strong absorption (Wo = 
2.5A) at this projected distance is very low. Further observations 
are needed in order to reveal the identity of the absorber towards 
PKS 0735+178. 
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